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Abstract 


This  program  outlines  a  three  step  process  to  develop 
nonlinear  composite  material  models,  and  incorporate  those 
models  in  a  general,  large  scale  structural  finite  element 
analysis  code.  The  steps  include  detailed  finite  element 
analyses  at  the  micromechanical  (fiber-matrix)  level  to  study 
material  response  and  nonlinear  mechanisms,  development  of  a 
simplified  material  model  that  is  representative  of  the 
behavior  observed  in  the  micromechanical  analyses,  and  the 
inclusion  of  the  simplified  model  in  the  large  scale  finite 
element  analysis.  Only  the  first  two  steps  are  addressed  in 
this  paper.  Micromechanical  finite  element  analyses  of 
AS4/3501-6  and  S2/3501-6  determined  stress-strain  responses 
for  in-plane  and  transverse  shear,  transverse  tension  and 
transverse  compression.  Variation  of  the  predicted  behavior 
with  failure  criteria  and  comparisons  to  experimental  data 
were  also  sho%m.  The  performance  of  a  simplified  material 
model  based  on  a  unit  cell  analogy  was  compared  to  the 
behavior  observed  in  the  micromechanical  analyses. 


Administrative  Information 

This  program  was  supported  by  the  DTRC  IR/IED  Program 
Office,  sponsored  by  ONR,  and  administered  by  Bruce  Douglas, 
DTRC  0112,  under  Work  Units  1-2844-220  and  1-2844-240. 

ZMTRODUCTIOK 

Composite  materials  can  exhibit  a  significant  nonlinear 
stress-strain  response  due  to  a  variety  of  mechanisms 
including  material  nonlinearities,  deunage,  and  interfacial 
debonding  [1].  These  nonlinearities  must  be  considered  for 
accurate  prediction  of  strength  or  stability  based  failure  of 
thick  composite  structures.  Through  thickness  stresses  can 
not  be  ignored  in  these  structures,  and  significant  material 
nonlinearities  are  expected  for  these  matrix  dominated  stress 
states . 

It  is  suspected  that  linear  material  assvimptions  have 
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several  effects  on  the  results  of  structural  analyses. 
Inaccurate  predictions  of  stresses,  strains,  deflections  or 
buckling  may  occur  if  material  softening  is  neglected.  These 
effects  are  often  counter  intuitive  with  composite  materials. 
On  the  one  hand,  a  nonlinear  structural  analysis  may  result 
in  more  critical  stresses,  strains  or  deflections  than  a 
linear  analysis.  On  the  other  hand,  stress  singularities 
predicted  with  linear  material  properties  may  be  relieved 
when  the  nonlinear  material  behavior  is  included  in  the 
structural  analysis  [2,3]. 

This  report  dociiments  the  development  of  a  three- 
dimensional  (3-D)  nonlinear  constitutive  modeling  approach 
for  composite  materials.  The  nonlinear  material  model  is 
based  on  a  combination  of  experimental  results  and 
micromechanical  finite  element  analysis  results. 

Incorporation  of  the  material  model  into  a  structural 
analysis  is  the  result  of  a  three  step  process: 

1)  development  of  the  detailed  3-D  nonlinear  response 
using  micromechanical  finite  element  analyses 

2)  formulation  of  a  simplified  material  model  whose 
behavior  matches  the  detailed  analysis  results 

3)  transition  of  the  model  to  a  general  purpose 
structural  analysis  finite  element  program. 

BACKGRODHD 

Early  efforts  in  composite  nonlinear  material  modeling 
focused  on  2-D  models  governed  by  continuum  mechanics 


formulat' jns.  Most  of  these  models  were  based  on  nonlinear 
elasticity  theory  [4,5,6].  other  approaches  included 
classical  incremental  plasticity  theory  [7]  and  endochronic 
plasticity  [8],  Griffin  et  al.  [3]  extended  the  effort  to 
three  dimensions  using  Hill's  anisotropic  formulation. 
Micromechanical  techniques  were  also  developed  to  address 
issues  ignored  by  continuum  mechanics  [9,10].  These  methods 
provided  detailed  information  on  matrix  and  fiber  stress 
distributions  for  any  combination  of  loading.  In  addition, 
they  could  be  used  to  evaluate  the  effect  of  constituent 
behavior  on  the  overall  composite  response. 

Although  recent  works  have  built  on  earlier  approaches  in 
classical  plasticity  [1]  and  endochronic  plasticity  [11], 
micromechanics  continues  to  be  a  strong  area  of  interest  due 
to  the  insights  it  offers  into  composite  material  behavior. 
Sun  and  Chen  [12],  Aboudi  [13]  and  Adams  [14]  have  developed 
unit  cell  models  of  varying  degrees  of  complexity  to 
represent  a  repeating  unit  within  a  unidirectional  lamina. 
These  models  are  then  used  to  predict  the  nonlinear  behavior 
of  a  laminate.  Sun  and  Chen  acknowledged  that  typical 
micromechanical  models  are  too  complex  to  use  in  the  analysis 
of  composite  structures.  They  used  results  (elastic-plastic 
stress-strain  curves  for  various  loading  combinations)  from 
their  detailed  micromechanical  model  to  develop  a  simpler, 
one  parameter  plasticity  approach  for  incorporation  in  a 
structural  analysis.  The  study  described  in  this  paper  uses 
a  similar  approach  where  a  simplified  unit  cell  model 
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developed  by  Pecknold  [15]  is  forced  to  match  the  behavior 
observed  in  detailed  micromechanical  finite  element  analyses. 


lATERIAL  MODEL  DEVELOPMEMT  STRATEGY 

Several  requirements  were  established  for  the  material 
model  in  this  progreun.  First,  the  model  must  be  suitable  for 
use  in  typical  composite  structural  analysis  tools.  This 
feature  provides  a  means  to  evaluate  material  constitutive 
effects  on  structural  performance.  However,  incorporation  of 
a  material  model  into  a  large  scale  structural  finite  clement 
program  limits  the  complexity  of  the  material  model  due  to 
the  large  number  of  times  the  material  model  is  invoked  in 
these  analyses  [12,15].  A  complex  material  model  would  result 
in  prohibitive  computer  run  times. 

Secondly,  the  nonlinear  behavior  must  be  characterized 
in  three  dimensions  for  any  given  loading  condition  by  the 
model.  Extension  beyond  the  2*'D  characteristics  of  classical 
lamination  theory  (CLT)  is  necessary  so  that  accurate  out-of¬ 
plane  material  behavior  can  be  included  in  structural 
analyses.  Characterization  of  the  overall  nonlinear  behavior 
covers  plasticity  and  damage  effects,  both  of  which  could 
have  significant  impact  on  structural  analysis  results 
[1,2,3]. 

Evaluation  of  the  stress-strain  response  for  any  given 
loading  condition  is  extremely  difficult  to  accomplish 
experimentally  due  to  the  inherent  problems  in  applying 
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combined  loads  on  material  test  specimens.  And  since  the 
material  behavior  is  nonlinear,  the  individual  uniaxial 
responses  cannot  be  superimposed.  Material  response  under 
multiaxial  load  conditions  is  particularly  useful  for 
structural  analyses  as  the  critical  areas  are  often  subject 
to  complex  stress  states.  A  model  that  can  provide  insight 
into  the  nonlinear  material  response  for  any  loading 
combination  can  be  used  in  lieu  of  testing. 

The  preceding  requirement  suggests  the  use  of  a 
micromechanical  material  model  where  the  constituent  (fiber 
and  matrix)  properties  can  be  related  to  the  composite 
properties.  Any  loading  combination  can  be  applied  to  the 
micromechanical  model  which  is  then  used  to  generate  the 
desired  composite  response  [12,13,15]. 

These  requirements  led  to  the  development  of  a  3-D 
nonlinear  material  modeling  scheme  that  is  conceptually 
similar  to  that  used  by  Sun  and  Chen  [12]  for  metal  matrix 
composites.  It  consists  of  three  steps:  development  of  the 
detailed  3-D  nonlinear  response,  formulation  of  a  simplified 
material  model,  and  transition  of  the  model  to  a  general 
purpose,  structural  analysis  finite  element  program,  such  as 
ABAQUS  [16],  through  a  user-written  material  (UHAT) 
subroutine . 

Detailed  Micromechanical  Analyses 

Micromechanical  methods  and  selective  experimentation 
are  used  to  develop  an  overall  understanding  of  the  3-D 
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material  Inelastic  response.  In  the  micromechanical 
approach,  composite  behavior  is  analyzed  using  some  type  of 
idealized  model  of  the  individual  fibers  and  surrounding 
matrix  material  [14].  The  idealized  model  may  be  as  simple 
as  a  Rule  of  Mixtures  (ROM)  formulation  [17]  or  as  complex  as 
a  finite  element  analysis  of  a  representative  volume  of  a 
unidirectional  composite.  The  latter  methodology  was  chosen 
for  this  program  as  it  most  closely  models  the  actual 
composite. 

wy02D,  a  2-D  finite  element  analysis  program  developed 
by  the  Composites  Materials  Research  Group  in  the  Department 
of  Mechanical  Engineering  at  the  University  of  Wyoming,  was 
selected  for  the  micromechanical  finite  element  anaylses. 

This  progreun  is  tailored  for  micromechanical  analyses  of 
continuous  fiber  reinforced  composite  materials  and 
incorporates  many  desirable  modeling  and  material  features. 
Since  the  source  code  is  provided,  the  program  can  be  easily 
modified  by  the  user. 

WY02D  represents  a  unidirectional  composite  as  a 
collection  of  unit  cells  consisting  of  a  fiber  embedded  in  a 
matrix  cube.  Using  a  generalized  plane  strain  approach,  this 
3-D  unit  cell  is  reduced  to  a  2-D  case.  Due  to  symmetry 
conditions,  only  one  quadrant  of  the  2-D  unit  cell 
representation  needs  to  be  modeled  (see  Figure  l) .  A  sample 
mesh,  including  fiber,  matrix  and  interface  elements,  is 
shown  in  Figure  2  [18].  Constant  strain  triangles  are  used 
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for  computational  speed  and  the  finite  element  grid  is 
automatically  scaled  to  reflect  changes  in  fiber  volume 
fraction. 

Although  a  rectangular  packing  array  is  assumed  for  the 
fibers,  any  fiber  packing  geometry  can  be  accomodated. 

Voids,  disbonds,  defects  and  degraded  interfaces  can  also  be 
modeled.  Using  a  technicnie  detailed  in  Reference  [14], 
displacement  boundary  conditions  are  applied  to  allow  the 
simultaneous  application  of  any  combination  of  loads. 
Temperature  and/or  moisture  content  changes  can  be  induced  to 
represent  such  phenomena  as  cooldown  induced  residual 
stresses,  moisture  induced  swelling  stresses,  and  moisture 
induced  constitutive  changes  in  the  matrix  material. 

Constituent  modeling  features  allow  the  use  of 
anisotropic  fiber  and  elastoplastic  matrix  properties.  The 
plastic  response  of  the  matrix  is  governed  by  the  Prandtl- 
Reuss  flow  rule  which  assumes  that  plastic  strain  is 
proportional  to  the  deviatoric  stress  [14,19].  The 
constitutive  equations  developed  with  this  flow  rule  require 
expressions  for  the  octahedral  shear  stress>strain  behavior. 
These  expresssions  are  obtained  from  a  Richard-Blacklock 
curve  fit  of  the  matrix  shear  stress-strain  data  as  shown  in 
Equation  1. 

T  »  (G€)/[1  +  (G€/To)"]^/*»  (1) 

r  «  shear  stress 
€  ^  shear  strain 

G  «  initial  tangent  shear  modulus 
TQ,n  *  curve  fitting  parameters 
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Various  failure  criteria  including  the  octahedral  shear 
stress,  maximum  normal  and  maximum  shear  [20]  can  be  chosen 
by  the  user  in  WY02D  to  determine  element  failures.  If  the 
element  stress (es)  exceeds  the  specified  criteria,  the 
element  stiffness  is  reduced  to  near  zero  and  the  strain 
energy  of  the  failed  element  is  redistributed  through  nodal 
forces . 

Small  modifications  to  the  program  were  made  by  the 
authors.  An  option  to  ignore  any  failure  criteria  in  the 
analysis  was  included  to  look  solely  at  matrix  plasticity 
effects  on  the  overall  nonlinear  stress-strain  response.  In 
addition,  the  interface  material's  tensile  strength  was 
set  equivalent  to  93%  of  the  matrix  dry  strength,  which  is 
halfway  between  the  matrix's  room  temperature  dry  and  wet 
strengths.  This  Interfacial  strength  was  chosen  based  on 
scanning  electron  microscopy  (SEM)  observations  of 
interfacial  failures  in  wet  and  dry  3501-6  composites.  The 
fiber  and  matrix  material  properties  used  for  analysis 
purposes  are  listed  in  Tables  1  and  2. 

Material  Model  ISoring  Model) 

Within  a  larger  scale  finite  element  analysis,  material 
model  calculations  must  be  performed  at  each  integration 
point  for  all  the  elements.  This  occurs  for  every  iteration 
(to  convergence)  within  a  load  increment  over  the  total 
nvimber  of  load  increments.  A  material  model,  such  as  WY02D, 
that  is,  in  itself,  a  detailed  finite  element  analysis  would 
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make  the  process  computationally  infeasible.  Therefore,  a 
technique  was  required  to  siimmarize  the  detailed 
micromechanical  results  in  a  simplified  form  so  that 
reasonedsle  computation  times  were  achieved. 

This  was  accomplished  using  a  material  model,  outlined  by 
Pecknold  [15],  whose  behavior  was  forced  to  resemble  that  of 
the  micromechanical  analyses  performed  with  WY02D.  In 
Pecknold 's  approach,  a  simplified  unit  cell  and  stress 
averaging  techniques  are  used  to  develop  the  response  of  a 
unidirectional  leunina.  The  spring  analogies  to  the  unit  cell 
configuration  in  Figure  3  are  based  on  the  variations  in 
load  path  for  in-plane  longitudinal  loads  and  shear/out-of¬ 
plane  loads.  Average  stresses  and  strains  within  the  unit 
cell  are  related  by  the  tangent  stiffness  for  the  unit  cell. 
The  tangent  stiffness  is  a  combination  of  the  tangent 
stiffness  parameters  of  the  fiber  and  resin  components  of  the 
iinlt  cell. 

As  in  WY02D,  the  fiber  is  considered  transversely 
isotropic.  In  contrast  to  the  elastoplastic  matrix  material 
model  in  HY02D,  the  matrix  shear  behavior  is  nonlinear 
elastic  and  expressed  in  a  Ramberg-Osgood  [21]  form. 

€  *  T/G  +  /5(Tq/G)  (T/Tq)**  (2) 

e  «  shear  strain 
r  «  shear  stress 

G  B  initial  tangent  shear  modulus 
«  curve  fitting  parameters 
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The  material  constants  used  in  the  spring  model  are  listed 
in  Tedsles  1  for  the  fiber  (same  as  WY02D)  and  Table  3  for  the 
3501-6  epoxy  matrix. 

Unidirectional  lamina  within  the  spring  model  are 
assembled  into  sublaminates,  the  smallest  repeating  unit  of 
the  overall  leuninate.  The  conditions  governing  the  assembly 
are  based  upon  classical  laminate  theory  with  one  extension. 
The  out-of-plane  stresses  at  each  leunina  interface  must  be 
continuous  across  the  interface  to  satisfy  equilibrium 
[22,23].  Development  of  the  sublaminate  tangent  stiffness  is 
dependent  upon  smearing  techniques  where  the  inhomogeneous 
l2uninate  is  replaced  by  a  homogeneous  anisotropic  material. 

Finite  Element  Computations 

Within  a  large  finite  element  program,  a  material 
svibroutine  is  typically  supplied  with  vectors  containing 
total  stresses  and  strains  at  the  beginning  of  the  load 
increment,  and  strain  increments.  Using  this  information, 
incremental  stress  updates  must  be  determined  within  the 
subroutine.  The  algorithm  chosen  for  these  updates  should  be 
computationally  efficient  due  to  the  large  number  of  times  it 
is  executed  in  the  analysis. 

This  study  used  an  algorithm  written  by  Hajali  [25] .  The 
algorithm  increments  strains  directly  by  taking  the  strain  at 
the  previous  increment  and  adding  on  the  strain  increment. 
Both  the  strain  at  the  previous  increment  and  the  strain 
increment  are  provided  as  input  to  the  material  sxibroutine 
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from  the  main  finite  element  program.  Micro  level  stress  and 
strain  increments  are  approximated  using  global  stress  and 
strain  increments  and  tangent  material  properties.  The  micro 
level  stress  and  strain  increments  are  used  to  calculate  new 
strain  totals  which  must  be  adjusted  in  an  iterative  process 
(modified  Newton-Raphson  method)  to  satisfy  equilibrium, 
compated^ility  and  nonlinear  constituent  equations.  Future 
refinements  of  the  stress  update  procedure  are  expected  to 
improve  computational  speed  and  convergence  time. 

RESULTS  RMD  DISCUSSION 

Detailed  Micromechanics  Model  Results;  Uniaxial  Loading 
Four  uniaxial  load  cases  were  studied  using  WY02D  for 
unidirectional  AS4/3501-6  and  S2/3501-6  laminates.  These 
included  in-plane  shear,  transverse  shear  (2-3  plane  shear) , 
transverse  tension  and  transverse  compression  as  shown  in 
Figure  4.  For  each  load  case,  the  stress-strain  behavior  to 
failure  was  predicted  using  different  failure  criteria 
(octahedral,  maximum  normal  stress,  and  maximum  shear 
stress) .  These  behaviors  were  compared  to  the  response 
predicted  with  no  failure  criteria,  i.e.  assuming  material 
plasticity  effects  only.  Figures  5  through  12  show  all  the 
generated  stress-strain  responses. 

When  availeUsle,  test  data  was  superimposed  upon  the 
predicted  behavior.  Figures  13  through  16  include  in-plane 
and  transverse  shear  data  from  the  losipescu  V-notched  beam 
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specimen  in  the  modified  Wyoming  test  fixture  [26]  and/or 
+45°  tension  [27]  test  for  unidirectional  AS4/3501-6  and 
S2/3501-6  l2uninates.  Each  of  these  methods  is  commonly  used 
to  determine  the  shear  modulus  in  the  linear  range.  The 
losipescu  specimen  is  also  used  for  shear  strength 
determination.  The  measured  nonlinear  response  is  quite 
different  for  each  specimen. 

No  attempt  was  made  to  force  the  predicted  behavior  to 
match  either  set  of  data  as  has  been  commonly  done  in  other 
micromechanical  approaches  to  account  for  uncertainties  in 
the  parameters  describing  the  fiber  and  matrix.  The 
predicted  response  is  based  on  a  unit  cell  subjected  to  a 
pure  and  uniform  shear  load.  In  reality,  the  losipescu  and 
+45°  tension  specimens  have  a  more  complex  stress  state  due 
to  geometry,  loading  effects,  and  damage  [28].  The  data  and 
predictions  may  not  match  depending  on  the  severity  of  these 
effects.  If  the  predicted  response  is  considered 
representative  of  the  true  material  response,  then  structural 
analyses  of  the  shear  specimens  using  the  nonlinear  material 
model  should  match  the  data  in  the  region  where  the  strain 
gages  were  mounted.  Of  course,  this  is  contingent  on 
accurate  loadings  and  boundary  conditions  in  the  structural 
analysis.  The  correlation  of  these  analyses  with  specimen 
data  will  be  documented  in  two  ensuing  reports. 

Transverse  tension  data  was  obtained  at  DTRC  using  an 
IITRI  specimen  [29].  Data  is  availed)le  for  S2/3501-6  only. 
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As  seen  in  Figure  17,  measured  and  predicted  initial  tangent 
moduli  correlate  well,  but  the  data  shows  more  softening  than 
predicted. 

Figures  18  and  19  show  transverse  compression  data  from 
Reference  30  for  AS4/3501-6  and  82/3501-6.  The  predicted 
octzdiedral  response  correlates  well  with  the  AS4/3501-6  data 
until  just  prior  to  failure.  The  transverse  compression  data 
for  S2/3501-6  is  highly  nonlinear,  more  so  than  any  of  the 
predicted  responses. 

In  the  preceding  comparisons  of  predicted  responses  and 
experimental  data,  it  appears  as  if  the  choice  of  failure 
criteria  is  dependent  on  the  particular  load  case.  However, 
these  comparisons  must  be  interpreted  carefully.  Direct 
comparisons  are  valid  only  to  the  extent  that  the  stress 
state  in  the  specimen  gage  section  is  uniaxial  and  uniform. 
Stress  analyses  of  the  specimens  incorporating  both  material 
nonlinearities  and  damage  induced  nonlinearities  would  yield 
more  accurate  comparisons. 

Micromechanics  Model  Results;  Multiaxial  Loading 

Any  combination  of  loads  can  be  applied  to  a  finite 
element  unit  cell  model  to  assess  the  effect  of  multiaxial 
loading  on  strength  and/or  stress-strain  behavior.  The 
difficulty  comes  in  obtaining  multiaxial  test  data  for 
verification.  One  comparison  of  data  and  predicted  response 
was  made  using  experimental  data  obtained  by  Davis  [31]  on 
boron/epoxy  tubes  subjected  to  combined  shear  (torque)  and 
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axial  compression.  Shear  stress-strain  curves  were  with  the 
tubes  under  three  levels  of  axial  compression  (  0,  1048,  and 
1572  MPa) .  Theoretical  simulations  of  this  load  state  were 
done  by  applying  in-plane  shear,  transverse  shear  and  axial 
compression  loads  to  the  unit  cell.  The  transverse  shear 
loads  represented  a  component  of  the  axial  compression  load 
on  the  tiibe  that  is  generated  due  to  imperfections  such  as 
fiber  waviness.  The  magnitudes  of  the  transverse  shear  loads 
were  calculated  numerically  by  Davis.  Figure  20  shows  the 
predicted  and  experimental  stress-strain  curves. 

Simplified  Material  Model  Rgsuj^g;  Uniaxial 
Loading 

Initially,  the  Ramberg-Osgood  parameters  describing 
the  matrix  in  the  spring  model  were  chosen  to  exactly  match 
the  Richard-Blacklock  model  of  the  matrix  shear  behavior  in 
HY02D.  The  predicted  matrix  shear  response  for  each  model  is 
shown  in  Figure  21.  The  in-plane  shear  response  for  a 
unidirectional  lamina  using  the  spring  model  was  then 
compared  to  that  generated  with  the  detailed  finite  element 
analysis  (WY02D) .  The  spring  model  approach  was 
significantly  stiffer  as  shown  in  Figure  22.  The  parameters 
describing  the  matrix  within  this  unit  cell  were  then 
adjusted  so  that  the  predicted  in-plane  shear  response 
matched  the  predicted  in-plane  shear  response  of  the  detailed 
FEA  performed  with  WY02D  for  both  AS4/3501-6  and 
S2glass/3501-6.  Figure  23  compares  the  shear  behavior  of  the 
matrix  constructed  with  these  parameters  to  the  matrix  in 
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Wy02D. 


The  best  correlation  between  the  spring  model  and  FEA 
results  for  in-plane  shear  was  obtained  with  P=l,  Tq-106  MPa, 
M-5  and  G»1.886  GPa  as  shown  in  Figures  24  and  25  for 
graphite  and  glass/ epoxy.  Using  the  same  parameters, 
comparisons  of  the  transverse  shear,  transverse  tension  and 
transverse  compression  responses  (unidirectional  leunina) 
using  the  spring  model  and  WY02D  are  shown  in  Figures  26 
through  31. 

The  AS4/3501-6  behavior  predicted  by  the  spring  model 
correlates  well  with  the  micromechanical  results  for  each 
load  case  without  a  change  in  spring  model  parameters  from 
load  case  to  load  case.  This  is  not  observed  for  the 
S2/3501-6.  The  parameters  that  are  used  to  correlate  the  in¬ 
plane  shear  behavior  do  not  result  in  matched  behavior  for 
the  remaining  load  cases.  Ho  attempt  was  made  to  address 
this  discrepancy  for  the  S2/3 501-6. 

Finite  Element  Computations 

Structural  analyses  performed  with  ABAQUS  using  the 
described  nonlinear  material  models  will  be  documented  in 
later  reports.  The  cases  studied  include  in-plane  shear 
specimens,  transverse  shear  specimens  and  composite 
cylinders. 
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CONCLUSIONS 


An  approach  has  been  implemented  to  incorporate  3-D 
nonlinear  material  behavior  into  structural  analyses.  The 
three  steps  of  this  approach  include  development  of  the 
detailed  3-D  nonlinear  response  using  micromechanical  finite 
element  analyses,  formulation  of  a  simplified  material  model, 
and  transition  of  the  model  to  general  purpose  structural 
analysis  finite  element  progreuns.  These  steps  provide  a 
range  of  methods  to  evaluate  nonlinear  behavior  from  the 
study  of  nonlinear  mechanisms  at  the  fiber-matrix  level  to 
the  assessment  of  nonlinear  material  effects  on  structural 
performance . 

Detailed  finite  element  micromechanical  analyses  form  the 
basis  for  the  constitutive  model  in  this  program.  They  offer 
a  \inigue  capability  to  evaluate  response  to  combined  loadings 
and  to  incorporate  the  use  of  failure  theories  to  evaluate 
the  effect  of  damage  on  composite  material  nonlinearity. 
Comparison  of  analysis  results  with  available  uniaxial  data 
showed  good  correlation  for  AS4/3501-6  laminates  loaded  in 
transverse  compression,  but  not  for  S2glass/3501-6  laminates. 
However,  transverse  tension  data  for  S2glass/3501-6  compared 
favorably  with  the  predicted  behavior.  Correlation  of 
analysis  with  in-plane  and  transverse  shear  data  will  be  done 
in  ensuing  reports  where  specimen  structural  analyses  are 
performed  with  the  nonlinear  material  models  developed  here. 
The  results  of  the  micromechanical  analyses  match  limited 
experimental  data  for  combined  compression  and  shear  loading. 
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Rigorous  finite  element  analyses  at  the  fiber-matrix 
level  are  too  complex  to  interface  directly  to  structural 
finite  element  analysis  codes.  A  simplified  material 
approach  called  the  spring  model  was  evaluated  for  this 
purpose.  The  performance  of  the  spring  model  correlated  well 
with  AS4/3501-6  micromechanical  analysis  results,  but  less  so 
with  the  S2glass/3501-6  results. 

A  subroutine  incorporating  the  spring  model  is  available 
for  ABAQUS  and/or  similar  finite  element  analysis  programs. 
Optimization  of  the  stress  updating  procedure  within  the 
subroutine  may  be  necessary  to  increase  computational 
efficiency. 

rUTURB  WORK 

It  was  desirable  to  incorporate  this  constitutive  model 
in  a  larger  scale  finite  element  program  as  quickly  as 
possible  to  determine  the  material  nonlinearities  critical  to 
structural  performance.  Based  on  this  decision,  WY02D  and 
Pecknold's  spring  model  were  chosen  as  they  were  readily 
available,  flexible  tools,  other  options  included  using 
comparable  micro-macro  material  models  or  writing  our  own. 

The  use  of  currently  available  tools  allowed  us  to  focus  on 
critical  unknovms.  There  are  several  areas  where  additional 
work  is  expected.  The  first  is  in  micromechanics.  It  is 
likely  that  the  failure  criteria  in  WY02D  will  be  replaced  by 
criteria  more  suitable  to  the  particular  matrix  system.  For 
example,  a  modified  Mohr's  failure  criteria  [20]  for  brittle 
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materials  with  different  tensile  and  compressive  strengths 
may  be  ideal  for  modeling  the  3501-6  epoxy  system.  This 
failure  criteria  is  more  complex,  in  itself,  requiring  a 
numerical  solution.  In  addition,  the  plasticity  model  may  be 
changed  to  reflect  new  materials. 

Other  ongoing  programs  are  generating  experimental  data 
on  strength  and  failure  mode  variations  in  wet  and  dry 
composite  materials.  The  insight  provided  in  these  programs 
on  matrix  and  interfacial  influences  will  be  incorporated 
into  the  detailed  micromechanical  analyses. 

Another  area  of  concern  regards  the  simplified  unit  cell 
model.  It  is  possible  that  this  model  may  not  track  some 
critical  behavior,  particularly  the  response  to  multiaxial 
loadings.  Methods  of  incorporating  failure  or  damage  into 
this  model  must  also  be  considered. 
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PROPERTY 

_ FIBER _ 

AS4 

S2GLASS 

Ei  (GPa) 

224 

88 

E2  (GPa) 

14 

88 

E3  (GPa) 

14 

88 

vi2 

.2 

.22 

^13 

.2 

.22 

''23 

.25 

.22 

Gi2  (GP®) 

14 

35 

Gi3(GPa) 

14 

35 

G23(GPa) 

5.6 

35 

Subscripts  indicate  direction 
1  Along  fiber 
2,3  Transverse  to  fiber 


T«bl«  1.  JIS4  grmpbltm  «nd  S2  glmm»  tibrnr  propertlM 
OMd  In  analysis 
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Property 

Epoxy  Resin 

3501-6 

Room  Temperature,  0%  Moisture 

E  (GPa) 

4.3 

G  (GPa) 

1.9 

V 

.34 

n 

1.604 

xq  (MPa) 

160 

Richard-Blacklock  Shear  Model 


Gy 

t  =  [  1  + 

X  Shear  stress 

Y  Shear  strain 

G  Initial  tangent  shear  modulus 

n,XQ  Curvature  parameters 


Tabl*  2.  3501-6  mpoxy  propaxtl**  ua*d  In  flnltn 

nil nnt  (inr02D)  ■icro— chani  cal  aMdallng 
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Property 

Epoxy  Resin 

3501-6 

Room  Temperature,  0%  Moisture 

K  (GPa) 

2.0 

G  (GPa) 

1.9 

P 

1.0 

N 

xq  (MPa) 

106 

Bulk  Modulus 
K=  E/(9-3E/G) 

Ramberg  Osgood  Shear  Model 
y  =  t/G  + 

T  Shear  stress 

y  Shear  strain 

G  Initial  tangent  shear  modulus 

P.N.tq  Curvature  parameters 


Tabl*  3.  3501-6  apoxy  propaztlaa  osad  in 

•priag  aodal 
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Finite  element  mesh  of  WY02D  unit  cell 
quadrant  (after  reference  18) 


3 


- 

1 

7  =  Fiber  Direction 


Unit  Cell  (V  =  Fiber  Volume  Fraction) 


For  component  1 1  f^or  components 

22.  12,  33.  23.  13 


Figure  3.  Spring  model  unit  cell  analogy  showing 
load  path  variations  between  fiber  and 
matrix  dominant  components 
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3 


In-Plane  Shear 
(1-2  Plane) 


Transverse 
Tension/Compression 
(1-2  Plane) 


Transverse  Shear 
(1-3, 2-3  Planes) 


1 


Directions 

1  Along  fiber 

2,3  Transverse  to  fibers 


fiqnxm  4.  UlcroMchanlcal  aodel  loadlxig  directions 
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IN-PLANE  SHEAR  STRESS  (MPA)  E  IN-PLANE  SHEAR  STRESS  (MPA) 


rigiix* 


WY02D  PREDICTED  IN-PLANE  SHEAR  BEHAVIOR 
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I  5.  iryo2D  px«dlet*d  ln-jplan«  abaar  bahavlor  for  AS4/3501-6 
nalag  no  failure  criteria  (material  plaatlclty  effects 
only),  and  the  octahedral,  normal  stresa,  and 

marimnm  shear  stress  failure  criteria 


WY02D  PREDICTED  IN-PLANE  SHEAR  BEHAVIOR 


IN-PLANE  SHEAR  STRAIN 


C.  1fy02D  predicted  In-plane  shear  behavior  for 

82glass/3501-6  using  no  failure  criteria  (material 

plasticity  effects  only),  and  the  octahedral, 

normal  stress,  and  marl mum  shear  stress  failure  criteria 


WY02D  PREDICTED  TRANSVERSE  SHEAR  BEHAVIOR 
AS4/3501-6 


rigur*  7.  inr02D  pradlctad  trmnmvmr—  mhmmx  bahavlor  for  AS4/3501-6 
ualag  no  failuro  crltorln  (antorlaa  plasticity  offoets 
only),  and  tha  oetahadral,  aariwai  noxaal  stxass,  and 
shaar  stzasa  fallnxa  erltarla 


WY02D  PREDICTED  TRANSVERSE  SHEAR  BEHAVIOR 
S2/3501-6 


Flgrura  8.  inf02D  pradletad  transvarsa  shaar  bahavlor  for 

S29lass/3501-6  using  no  fallnra  erltarla  (utarlal 
plasticity  affacts  only),  and  tha  oetahadral,  asTlimui 
nozaal  strass,  and  ■si’lna  shaar  strass  fallura  erltarla 
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TRANSVERSE  STRESS  (MPA)  8  TRANSVERSE  STRESS  (MPA) 


WY02D  PREDICTED  TRANSVERSE  TENSION  BEHAVIOR 
AS4/3501-6 
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9.  Wyo2D  predicted  transv<«rstt  tension  behavior  for  hS4/3501-€ 
using  no  failure  criteria  (aaterial  plasticity  effects 
only),  and  the  octahedral,  nomal  stress, 

■arlue  shear  stress  failure  criteria 


WY02D  PREDICTED  TRANSVERSE  TENSION  BEHAVIOR 
S2/3501-6 


Figure  10.  irr02D  predicted  transverse  tension  behavior  for 

S2glass/3501-6  using  no  failure  criteria  (siaterlal 
plasticity  affects  only) ,  and  the  octcthedral,  siaxiauiD 
norswl  stress,  and  ■axleum  shear  stress  failure  criteria 
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Flgor*  11.  I(Y02D  pr«dlet«d  trans^rs*  coaiprasBlon  bahaTlor  for 
AS4/3501-6  ualag  no  fnlluro  erltorla  (natarlal 
plasticity  of fact a  only),  and  tho  octahodral, 

narlow  noroal  stross,  and  ■arlMi  aboar  stross 

fallozo  erltorla 


TRANSVERSE  COMPRESSIVE  STRAIN 


riguro  12.  1(y02D  prodlctod  transvorso  coaiprosslon  boba^lor  for 

S2glass/3501-6  using  no  falluro  erltorla  (aatorlal 
plasticity  offsets  only) ,  and  tbo  octahodral, 

— noxaal  stross,  and  ■a¥^^m^,  ahoar  stross 
falluro  erltorla 
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IN-PUNE  SHEAR  STRESS  (MPA)  H  IN-PUNE  SHEAR  STRESS  (MPA) 


IN-PLANE  SHEAR  BEHAVIOR 
PREDICTED  VS  DATA 
AS4/3501-6 
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I  13.  Conparison  of  prodlctod  AS4/3501-6  In-plaa*  shoar 

bahavlor  with  Zoalpaaeo  and  ±45^  tansloa  shaar  data 
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F:  /urm  14.  Coaparlson  of  pxadletod  82glasa/3501-6  In-plana  shaar 
bahawior  with  loslpasco  aad  ±45°  tanslon  shaar  data 


IN-PUNE  SHEAR  BEHAVIOR 
PREDICTED  VS  DATA 
S2/3501-6 
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TRANSVERSE  SHEAR  STRESS  (MPA)  3  TRANSVERSE  SHEAR  STRESS  (MPA) 


TRANSVERSE  SHEAR  BEHAVIOR 
PREDICTED  VS  DATA 
AS4/3501-6 


15.  Coaparlson  ot  pradletad  AS4/3501-6  tTmnmvrmm  mhmmr 
b«ha^lor  with  Zoslpwseu  shwar  data 


TRANSVERSE  SHEAR  BEHAVIOR 
PREDICTED  VS  DATA 


rigura  1€.  Coa^parlson  of  pradletad  S2glas8/3501-6  traaswarsa 
•haar  bahavlor  with  loslpaaeu  ahaar  data 
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TRANSVERSE  TENSION  BEHAVIOR 
PREDICTED  VS  DATA 
S2/3501-6 


Flgor*  17.  Coap«rl«on  of  prodlctod  S2gXass/3501-6  tranawso 
tonalon  b«h«Tlor  with  IZTRI  data 
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TRANSVERSE  COMPRESSION  BEHAVIOR 
PREDICTED  VS  DATA 


Flgur*  18.  Coaiparlson  of  prodicted  Jk84/3501-6  tranarorso 
eoivrasslon  bahavlor  with  data 


TRANSVERSE  COMPRESSION  BEHAVIOR 
PREDICTED  VS  DATA 
S2/3501-6 


riguro  19. 


Coaiparlson  of  prodlctad  S2glaaa/3501-6  transvaraa 
eoaipsosslon  baharlor  with  data 
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SHEAR  STRESS  (MPA 


COMBINED  COMPRESSION  AND  TORSION  LOADING 
PREDICTED  VS  DATA 


SHEAR  STRAIN 


rigur*  20.  Coapulson  of  prodlctod  boroa/«poxy 

coaprosslon/abaar  loading  bahavlor  with  data  from 
Safaranoa  31 
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MATRIX  SHEAR  BEHAVIOR 
3501-6 


Figure  21.  Conparlson  of  matrix  shear  behavior  using  08g>lrlcally 
based  swdel  In  inr02D,  and  the  spring  model  with  matrix 
P»rmmmtmTm  chosen  so  as  to  match  the  m02D  model 


PREDICTED  IN-PLANE  SHEAR  BEHAVIOR 
SPRING  MODEL  VS  MICROMECHANICS 


Figure  22.  Comparison  of  predicted  AS4/3501-6  In-plane  shear 
behavior  using  matching  matrix  models  la  irr02d 
and  the  spring  model 
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SHEAR  STRESS  (MPA) 


MATRIX  SHEAR  BEHAVIOR 
3501-6 


SHEAR  STRAIN  (M/M) 


rigon  23. 


of  Mtrlz  shoar  bohaTlor  using 
siiplrlcslly  bsssd  aodsl  in  inr02D  and  xarlsad 
■atrlx  paraaatars  in  tha  spring  nodal 


39 


IN-PLANE  SHEAR  STRESS  (MPA)  f!  IN-PLANE  SHEAR  STRESS  (MPA) 


PREDICTED  IN-PLANE  SHEAR  BEHAVIOR 
MICROMECHANICS  VS  SPRING  MODEL 
AS4/3501-6 
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IN-PLANE  SHEAR  STRAIN 


gur*  24. 


1fy02D  and  spring  aodal  pradletad  ln-plnn« 
■hnnr  bahnvlor  tor  AS4/3501-6 
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TRANSVERSE  SHEAR  STRESS  (MPA)  TRANSVERSE  SHEAR  STRESS  (MPA) 


PREDICTED  TRANSVERSE  SHEAR  BEHAVIOR 
MICROMECHANICS  VS  SPRING  MODEL 
AS4/3501-6 


'Igure  26.  irr02D  and  apxlxtg  modal  pxmdlctad  tranwarsa 
shaar  bahavlor  tor  AS4/3501-6 


PREDICTED  TRANSVERSE  SHEAR  BEHAVIOR 
MICROMECHANICS  VS  SPRING  MODEL 
S2/3501-6 


Flgura  27 .  irf02D  and  aprlag  a»dal  pradietad  transvaraa 
•haar  baha^lor  tor  S29lasa/3501-6 
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TRANSVERSE  STRESS  (MPA) 


PREDICTED  TRANSVERSE  TENSION  BEHAVIOR 
MICROMECHANICS  VS  SPRING  MODEL 
AS4/3501-6 


Flguz#  28 ,  1fy02D  axul  spring  pcwUcfcAd  tnnsvsf 

tension  bahsTlor  for  JkS4/3501*S 


PREDICTED  TRANSVERSE  TENSION  BEHAVIOR 
MICROMECHANICS  VS  SPRING  MODEL 


Flgfuro  29.  irr02D  and  aprlng  nodal  pradletod  tranavarsa 
tanalon  bahavlor  for  82glass/3501-€ 
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PREDICTED  TRANSVERSE  COMPRESSION  BEHAVIOR 
MICROMECHANICS  VS  SPRING  MODEL 


rigor*  30.  1fy020  and  aprlxig  modal  pr*dlct*d  tran*o*rs* 

eoapr***lon  bahavlor  tor  AS4/3501-'6 


PREDICTED  TRANSVERSE  COMPRESSION  BEHAVIOR 
MICROMECHANICS  VS  SPRING  MODEL 


rigor*  31.  wy02D  and  aprlng  modal  pradlctad  traaaoara* 
eomprosslon  bahaTlor  for  S2glaaa/3501-6 
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